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Abstract 

New recursive equations designed for the G/G/m queue are pre- 
sented. Tliese equations describe the queue in terms of recursions for 
the arrival and departure times of customers, and involve only the op- 
erations of maximum, minimum and addition. 
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1 Introduction 

As a representation of dynamics of queueing systems, recursive equations 
have been introduced by Lindley in his classical investigation of the G/G/\ 
queue [T] . In the last few years the representations based on recursive equa- 
tions have been extended to a variety of systems which consist of single- 
server queues and operate according the first-come, first-served (FCFS) 
discipline. Specifically, there are the equations designed for open and close 
tandem queues with both infinite and finite buffers [21 [3l H]. Recursive 
equations have been also derived to represent more complicated systems of 
G/G/1 queues, including acyclic fork-join networks O |6] and closed net- 
works with a general deterministic routing mechanism [71 El |9] . 

Recursive equations find a wide application in recent works on both 
analytical study and simulation of queueing systems. As an analytical tool, 
they were exploited to investigate system performance measures [21 [71 [10] 
and estimates of their gradients [3 |S1 121 [HJ • In [5| [S] recursive equations 
based representations have provided the means for establishing the stability 
conditions and deriving bounds on system performance in a class of queueing 
systems. Finally, these representations made it possible to develop efficient 
algorithms of queueing systems simulation [3, ^ as well as powerful methods 
of estimating gradients of system performance measures [71 [9l [TT]. 
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In the queueing theory, the G/G/m queue is normally represented by 
using the recursions introduced by Kiefer and Wolfowitz in [12]. Expressed 
in rather general terms, these recursions may be inconvenient to use if an 
explicit form of representation is required. The purpose of this paper is 
to present new recursive equations designed for the G/G/m queue. These 
equations describe the queue in terms of the arrival and departure times 
of customers, and involve only the operations of maximum, minimum and 
addition. 

The rest of the paper is organized as follows. In Section 2 we briefly 
describe both the G/G/1 and G/G/m queues, and give preliminary analysis 
of their representations. Section 3 includes two technical lemmas which offer 
useful representations for the ordering operator. Finally, in Section 4 we 
present our main result providing an explicit recursive representation for 
the G/G/m queue. 

2 Definitions and Preliminary Analysis 

Recursive equations as a representation of queueing systems were originally 
expressed in terms of waiting times [TJ [T^] . Equations of this classical type 
continue in use (see, e.g., [TO]); however in many recent works [2 [3l IU 
[3 El El El El E] one can find the equations describing system dynamics 
through recursions for arrival and departure times. The last approach to 
the representation of queues may be considered as a useful generalization of 
the classical one, and it will be applied below to derive recursive equations 
for the G/G/m queue. 

We start with the equations for the G/G/1 queue, which provide the ba- 
sis for representing more complicated systems including the G/G/m queue. 
To set up these equations, consider a single server queue with infinite buffer 
capacity and the FCFS queue discipline. Denote the interarrival time be- 
tween the kih. customer and his predecessor by a^, and the service time 
of the feth customer by Tk- We assume > and > 0, for any 
k = 1,2,... Furthermore, let be the kth. arrival epoch to the queue, 
and Dk be the kih. departure epoch from the queue. As is customary, the 
sequences {ak}k>i and {Tk}k>i are assumed to be given, whereas {Ak}k>i 
and {-DA:}fc>i are considered unknown. Finally, provided that the queue 
starts operating at time zero, it is convenient to set = and Dk = for 
ah A; < 0. 

Using these notations, the recursive equations for the G/G/1 queue may 
be written as [21 El S] 

Ak = Ak^i + ak 

Dk = {Ak V Dk-i) + Tk, 

where V denotes the maximum operator, k = 1,2, . . . 
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The first equation in ([T]) is trivial. To understand the second recursion it 
is sufficient to see that the term AfcV-Dfc_i determines the service initiation 
time for the kth customer. Clearly, the service of this customer may be 
initiated not earlier than he arrives at the server. If the A:th customer 
finds the server busy, he has to wait until service of the (A; — l)st customer 
completes. In other words, the time of the kth initiation of service coincides 
with the maximum out of Af^ and -Dfc-i • 

Taking ([1]) as the starting point, we now turn to the discussion of multi- 
server queues. Consider a queueing system with an infinite buffer and m, 
1 < m < oo, servers operating in parallel. When a customer arrives, he 
occupies any one of those servers which are not busy. If there are no free 
servers at his arrival, the customer takes his place in the buffer and starts 
waiting to be served. His service begins as soon as any one of the servers is 
unoccupied, provided that all his predecessors have left the buffer. 

We may extend all the definitions introduced above to this queueing 
system. Note, however, that in the multi-server queue the kth departure 
time may not coincide with the completion time of the A;th customer. One 
can consider the first customer as an example. Being the first to initiate 
service, he may prove to be the kth to depart from the system, for any 
A: > 1 . To take account of the distinction between these times, we further 
introduce the notation Ck for the completion time of the kth customer, 
k = l,2,... 

Since upon their service completions the customers immediately leave 
the system, it is easy to see that both the sequences of completion times 
and departure times are constituted by the same elements. In fact, the 
sequence {Dk}k>i is simply the sequence {Ck}k>i arranged in ascending 
order. Therefore, an operator which produces ordered values is required to 
express the relation between these sequences. 

The idea to apply some ordering operator in recursive equations based 
representations of the G/G/m queue had its origin in [12j. However, this 
operator, as it has been introduced in [K] and occurs in recent works (see, 
e.g., [6]), is expressed in general terms rather than in an explicit form. It will 
be shown in the next section how the ordering operator may be represented 
in closed form function of values being ordered. 

3 Representations of the Ordering Operator 

Let {vj^^^i = {ri, . . . ,r„} be a finite set (sequence) of real numbers. Sup- 
pose that we arrange its elements in order of increase, and denote the kth 
smallest element by r"^^ . If there are elements of an equal value, we count 
them repeatedly in an arbitrary order. Finally, we introduce the notation A 
for the minimum operator, and set r^u\ = —oo for all k <0. 



3 



<^= A (^.iV-.-Vr,J. (2) 



Lemma 1. For each k = 1, . . . ,n , the value of r^.^ is given by 
(k) = 

l<ji<-<jfc<n 

The proof of this statement can be found in Now suppose that a 
new element r„+i is added to {rj}"^^. For the expanded set {rj}^^l, we 
denote the kth smahest element by i^"^^^ ■ The next lemma is intended to 
relate the ordered values of the set {rj}"^| to those of {rj}^^-^ . 



Lemma 2. For each k = 1, . . . ,n , it holds 

(k) 

(fc) 

(fc) 



r«+i=rf,)A(rf,_,)Vr„+i)- (3) 
Proof. To derive ([3]), we first apply Lemma 1 to represent r?j\^ in the form 



^-'= A (^ii V---Vr,J. 



l<ii<-<Jfc<n+l 

The terms being minimized may be rearranged to write this relation as 



^(fc) 



l<ji<-<jfc<n / \l<ji<-<jfe=n+l 



Clearly, the first term on the right-hand side of the above relation just 
represents r"^^ . It remains to factor out the common member rj^, = r„+i 
from the second term, and to apply Lemma 1 once again so as to get Q: 



^W'=^W^ A (^.-1 V--- Vr,,^_JVr„+i 

\i<ii<---<ifc-i<" 

= A(r[^fe_i) Vr„+i). □ 

To conclude this section, note that the above representations of the or- 
dering operator involve only the operations of maximum and minimum, and 
appear to be particularly suited for use in the recursive equations under 
discussion, which are actually expressed in terms of similar operations. 



4 Recursive Equations for the G/G/m Queue 

We are now in a position to present the main result of the paper. We start 
the derivation of the equations representing the G/G/m queue with the 
observation that the arrival process in this queue system is no different from 
that in the G/G/1. Therefore, the first equation at ([1]) 

Ak = Ak-i + ak 
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remains unchanged. 

Obviously, to calculate the completion time of the customer which is 
the kth to arrive to the system, one has to add his service time to 
the time of his service initiation. Similarly as in the above analysis of the 
G/G/1 queue, let us examine the possibilities for this customer to initiate 
his service. Firstly, the customer may find one or more servers free at his 
arrival. In this case he receives service immediately at the time Aj, . 

Suppose now that at the arrival of the fcth customer all the servers 
are found to be busy. If there are no other customers waiting for service, 
he occupies that server which becomes free earlier. Clearly, the customers 
being served unoccupy servers according to the sequence of the departure 
times Dk-m, ■ ■ ■ , Dk-i ■ We may therefore conclude that, as the earliest 
time in this sequence, Dk_rn just represents the service initiation time of 
the feth customer. It is easy to see that this conclusion remains the same 
in case there are previous customers still waiting for service, when the /cth 
customer arrives. Finally, for the completion time of this customer, we have 

Ck = {Ak V Dk-m) + Tk- 

As an important consequence of the above equation, one can state that 
Cfc > Dk^^. By renaming the indices, this relation may be rewritten as 

Dk < Ck+m- (4) 

To represent the G/G/m queue completely, we have to define the depar- 
ture time Dk by a suitable equation. One can conclude from the discussion 
in Section 2 that coincides with the klh. smallest element of the set 
{Cj}j>i. Moreover, as it follows from (jH), it is sufficient to examine only 
the finite subset {Ci, . . . , Ck+m-i} so as to determine D^. Let C^k)^~^ be 

the fcth smallest element of this set. We may now define D^. = C'(^^™'~^ • 

By applying ^ and (H]), and taking into account that Dfe-i = C^k-l)'^^ 
we successively get 

'-'{k) - ^(fc) ^^(k-l) ^ '^k+m-1) 

= C^i^"^ ^ A (Dfc-l V Ck+m-l) = C'^'fe)^™' ^ A Ck+m-1- 

Finally, using we arrive at 

Dk = CffcJ™-' = /\ (C,, V ... V Gj, ) A Gk+m-i- 

l<ji<---<jk<k+m-2 

We summarize the above results as follows. 
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Theorem 3. With the notations introduced in Section 2, the dynamics of 
the G/G/m queue is defined by the recursive equations 



Ak = Ak-i + Ok 

Ck = {Ak V Dk-m) + Tk 



(5) 



A 



(c,, V . . . V A a 



:+m— 1 ) 



l<ii<-<ifc<fc+m-2 



for each k = 1,2, . . . 

In conclusion, let us consider two consequences which follow immediately 
from the theorem. Clearly, if m = 1 , it results from the last equation at ^ 
that Ck = Dk, k = 1,2, . . . In this case, the above equations are reduced to 
the recursions ([1]) representing the G/G/1 queue. 

As another consequence of the obtained representation, one may produce 
the equations describing the G/G/2 queue. Putting m = 2 in ([5]), it is not 
difficult to arrive at 
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